Jet-induced gauge field instabilities in the quark-gluon plasma: A kinetic theory 
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We discuss tiie properties of the collective modes of a system composed by a thermalized quark- 
gluon plasma traversed by a relativistic jet of partons. The transport equations obeyed by the 
components of the plasma and of the jet are studied in the Vlasov approximation. Assuming that 
the partons in the jet can be described with a tsunami-like distribution function we derive the 
expressions of the dispersion law of the collective modes. Then the behavior of the unstable gauge 
modes of the system is analyzed for various values of the velocity of the jet, of the momentum of 
the collective modes and of the angle between these two quantities. We find that the most unstable 
modes are those with momentum orthogonal to the velocity of the jet and that these instabilities 
appear when the velocity of the jet is higher than a threshold value, which depends on the plasma 
and jet frequencies. The results obtained within the Vlasov approximation are compared with the 
corresponding results obtained using a chromohydrodynamical approach. The effect we discuss here 
suggests a possible collective mechanism for the description of the jet quenching phenomena in heavy 
ion collisions. 
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I ■ I. INTRODUCTION 

, One of the interesting aspects of ultrarelativistic heavy-ion coUisions is that the properties of the produced matter 
can be studied employing high px partons generated by hard scatterings in the initial stage of the collision. Such high 
Pt partons behave as hard probes unveiling some properties of the system and providing evidences of the production 
• of a thermalized quark-gluon plasma (QGP) m. When the jet of partons travels across the medium it loses energy and 
^C) ' degrades, mainly by radiative processes (see (2| for reviews). The energy and momentum of the jet are absorbed by 
0^ ■ the plasma and result in an increased production of soft hadrons in the direction of propagation of the partons. The 
OO ■ description of this process is currently undertaken employing various models, where perturbative QCD is supplemented 
i with medium- induced parton energy loss 0] or where the AdS/CFT correspondence [1] is employed. 
' In a previous paper [5| we have proposed a novel mechanism for describing how the jet loses energy and momentum 
' while traveling in a thermally equilibrated quark-gluon plasma. Since the jet of particles is not in thermal equilibrium 
with the QGP it perturbs and destabilizes the system inducing the generation of gauge fields. Some of these gauge 
modes arc unstable and grow exponentially fast in time absorbing the kinetic energy of the jet. A basic assumption 
in was that that both the plasma and the jet can be described using a fluid approach. This approach developed in 
@ has been derived from kinetic theory expanding the transport equations in moments of momenta and truncating 
^ , the expansion at the second moment level. The system of equations is then closed with an equation of state relating 
' pressure and energy density. The fluid approach has several advantages with respect to the underlying kinetic theory. 
The most remarkable one is that one has to deal with a set of equations much simpler than those of kinetic theory. 
Then one can easily generalize the fluid equations to deal with more complicated systems. This is a strategy that has 
been successfully followed in the study of different dynamical aspects of non-relativistic electromagnetic plasmas ■ 
In the present paper we consider the same setting of our previous paper Q , i.e. a static plasma traversed by a 
relativistic jet, but we use kinetic theory instead of the fluid approach. The appearance of filamentation instabilities 
due to hard jets was flrst pointed out in Ref. Q. Here we generalize that study, also with the aim of preparing the 
ground for the computation of the jet energy loss. However, in the present paper we content ourselves with the study 
of the collective modes of the system, postponing to a future publication [9| the computation of the jet energy loss. 

Transport theory provides a well controlled framework for studying the properties of the quark-gluon plasma in 
the weak coupling regime, 5 <C 1. Indeed it is well known that the physics of l ong distance scales in an equilibrated 
weakly coupled QGP can be described within semiclassical transport equations |lOl [ill. [T^. In this approach the hard 
modes, with typical energy scales of order T, are treated as (quasi-) particles which propagate in the background of 
the soft modes, whose energies are equal or less than gT, which are treated as classical gauge fields. This program 
has been very successful for understanding some dynamical aspects of the soft gauge fields in an almost equilibrated 
QGP [H, The same approach has been recently used to study the development of gauge field instabilities when 
the distribution function of the hard modes is far away from equilibrium, and it is anisotropic in momenta space (see 



2 



Ref. [l5| and references therein) . We will also use the semiclassical transport approach in order to study how the soft 
gauge field dynamics of a thermally equilibrated QGP is perturbed by a relativistic jet of particles. 

The hard modes of the QGP can be described with two different transport equations. One can treat the color 
charges as classical variables [13, HB] or one can treat color as a quantum degree of freedom [l^, . We will employ 
the latter formulation. However since we will study the gauge fields dynamics in the linear response approximation 
employing the former formulation we would have obtained the same results. Indeed in the leading response analysis 
the two possibilities lead to the same answers. We will also compare the results for the growing rates of the unstable 
modes obtained in transport theory with the analogous results we obtained with the fluid approach in Ref. 0. 

This paper is organized as follows. In Section we review the kinetic theory approach in the so-called Vlasov 
approximation and specialize the system of equations to the case of interest of a plasma traversed by a jet of particles. 
In Section [IIII we consider separately a thermally equilibrated plasma and a jet of particles with a so-called tsunami- 
like distribution function. In both cases we derive the dispersion laws of the collective modes and find that the two 
systems are stable. In Section IIVI we show that in a system composed by a thermalized plasma traversed by a jet 
of particles unstable collective modes appear, and we determine their growing rate. We discuss separately the cases 
where the momentum of the collective mode is collinear with the velocity of the jet, orthogonal or at an arbitrary 
angle. In Section |V] the results obtained within the kinetic theory approach arc compared with the results of the fluid 
approach. In Section IVII wc discuss our results and comment on how our study could be improved to make contact 
with heavy ion phenomenology. 



II. KINETIC THEORY FOR COLORED PARTICLES IN THE VLASOV APPROXIMATION 

In this Section we review the transport equations obeyed by the distribution functions of colored particles [13, [i3 
that represent the starting point of our analysis. 

The distribution functions of quantum colored particles are hermitian matrices, whose dimensionality depends on 
the color representation of the particle. They are not gauge invariant quantities, but rather transform covariantly 
under a gauge transformation. Quarks and antiquarks belong to the fundamental representation of SU (3)c and the 
corresponding distribution functions Q(j),x) and Q{p,x) are 3x3 matrices in color space. Gluons arc in the adjoint 
representation of SU (3)c and their distribution function G{p, x) is a 8 x 8 matrix in color space. Without loss of 
generality we will consider that the particles in the jet are quarks, and we will denote the corresponding distribution 
ftmction as Wjctip, x), separating explicitly that component from Q(p,x). The reason for doing so is that we will 
consider very different initial conditions for the thermally equilibrated quarks in the plasma and for the quarks in 
the jet. Also, one should note that the coUisional rates will probably be very different for these two systems of 
quasiparticles. 

The distribution function of quarks, antiquarks, gluons and of the particles of the jet satisfy the following transport 
equations 

p'^D^Q{p,x) + ^p^ {F^,{x),d;Q{p,x)} = C, (1) 
P^-D^Qip, - I {F^.{x), d;Q{p, x)} = C, (2) 
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pf'D^WietiP, + I p'' {F^,{x),d;Wi,tiP, x)} = Cw (4) 



P^V^G{p,x) + {T^,{x),d;G{p,x)} = Cg , (3) 



With {..., ...} we denote the anticommutator. dp is the four-momentum derivative and g is the QCD coupling constant. 
The covariant derivatives and act as 

Df^^df,- ig[Af,{x), ... ] , V^.^d^- ig[Af,{x), ... ] , 

with = A^{x)t'^ and = Ai;^{x)T°-, and and are SU{3) generators in the fundamental and adjoint 
representations, respectively. The strength tensor in the fundamental representation is i*),^ = d^j,Ay—duA^j^—ig[Af^, Ay\^ 
while T^iy denotes the field strength tensor in the adjoint representation. 

In Eqs. (HI),©, © and © C, C, Cg and Cw represent the colhsion terms. For time scales shorter than the mean free 
path time the collision terms can be neglected, as typically done in the so-called Vlasov approximation. The knowledge 
of the distribution function allows one to compute the associated color current, which in a self-consistent treatment 
enters as a source term in the Yang-Mills equation. Therefore in this approximation the different components of the 
system formed by the plasma and the jet interact with each other only through the generated average gauge fields. 
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Assuming that the quark-gluon plasma is initially in a colorless and thermally equilibrated state, we can describe 
the small deviation from equilibrium of the distribution functions of quarks, antiquarks and gluons as 

QiP, x) = fpoiPa) + ^QiP^ x) ' QiPi x) = fFDiPo) + ^Q{P, , G{P, = /b£(po) + SG{p, x) , (5) 

where 

are the Fermi-Dirac and Bose-Einstein equilibrium distribution functions. 

We will also study color fluctuations of the initial colorless distribution function of the jet. Thus 

M^jetb, X) = /jet(p) + <5M^jctb, X) . (7) 

For the initial jet distribution function we will consider a colorless tsunami-like form [Toj 

fj,t{p)=nu° S^'^{p~Au) , (8) 

that describes a system of particles of constant density n, all moving with the same velocity = (u°,u) = 7(1, v), 
where 7 is the Lorentz factor. The parameter A fixes the scale of the energy of the particles. 

Although this distribution function is adequate for describing a uniform and sufficiently dilute system of particles, 
it would be very interesting to extend our analysis to more complicated forms. Using more involved distribution 
functions on the one hand would lead to a more accurate description of the jets relevant for heavy ion phenomenology, 
when the density of particles composing the jet is not uniform and there is a spread in momentum. However, on the 
other hand this would complicate the dispersion laws of the collective modes. For this reason we leave such analysis 
to a future work. 

In the Vlasov approximation one can compute the contribution to the polarization tensor of the particles of species 
a (where a refers to quarks, antiquarks, gluons or the partons of the jet) as 

Kb,aik) = -9 CpSab J faip) ^^-^^^ , (9) 

where a, b are color indices and Cp is the value of the quadratic Casimir associated with the particle specie a which 
takes values 1/2 and 3 for the fundamental and adjoint representations, respectively. The momenta measure is defined 
as 



(2^ 



2e(po)<5(p^ ~ mi) , (10) 



where TOq, is the mass of the particle of specie a. For simplicity we assume that the particles belonging to the plasma 
are massless. Instead the particles of the jet have a non- vanishing mass. 

When fa{p) is a thermal equilibrated distribution function, Eq. (O reduces to the form of the hard thermal loop 
(HTL) polarization tensor [13 Ullj [13] ■ However for the tsunami-like distribution function, Eq. ([S]), the polarization 
tensor obviously takes a different form. 

The gauge fields obey the Yang-Mills equation 

D.F^-'ix) = SjUx) = Sj;{x) + Sj^,,{x) , (11) 

where we have defined 

Sj^ix) fp'' [SQiP^^) - SQ{p,x) + 2t-Tt[T-SG{p,x)]] , (12) 



which describes the plasma color current, and 



SjUx)^-^ I p''SWi.t{p,x), (13) 



which describes the fluctuations of the current associated with the jet. 
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Equation ([TT|) together with Eqs. (IT]|4]) form a set of equations that has to be solved self-consistently. Indeed the 
gauge fields which are solutions of the Yang-Mills equation enter into the transport equations of every particle species 
and, in turn, affect the evolution of the distribution functions. 

We will show that in the presence of a fast moving jet of particles some unstable gauge modes appear, as first pointed 
out in Ref. @, and also discussed using a chromohydrodynamical approach in our previous paper However it is 
well-known that if there is no jet of particles propagating in the plasma, then all the gauge collective modes are stable 
[20| . This is due to the fact that the HTL approximation is obtained employing a distribution function corresponding 
to a state of stable minimum for the plasma. In the next Section we will derive the equation of the collective modes 
of a thermally equilibrated QGP (where there is no jet of particles). Then we will derive the dispersion laws for the 
collective modes of a jet of particles with a tsunami-like distribution function and show that also in this case the 
collective modes are stable. 



III. COLLECTIVE MODES IN THE EQUILIBRATED QGP AND IN THE JET 



The collective modes of a thermally equilibrated plasma are obtained from the knowledge of the HTL polarization 
tensor. We briefly recall some results of the HTL approach and derive the dispersion laws that describe the evolution 
of gauge collective modes [20| . 

The dielectric tensor is defined as 



H'^ 



(14) 



where in order not to overcharge the notation we dropped the color indices. For the QGP in the HTL approximation 
the dielectric tensor is given by [lO] 



where 



eL(cc',k) = 1 + -^ 



eT(w,k) 



2k uo 



-r^eL(^,k) , 





k ^ 




i- 




uj + k 


\2k 





(15) 

(16) 
(17) 



are the longitudinal and transverse (with respect to k) components of the dielectric tensor of the plasma, composed 



by gluons and massless quark of Np flavors, and = ^g^T^{3 + Np/2) is the plasma frequency squared. 



The dispersion laws arc obtained by solving the equation 



det 



k^6'^ -k'k^ -oj^e'Jik) 



= 0. 



that leads to 



0, 
, 



(18) 



(19) 
(20) 



for the transverse and longitudinal modes respectively. In general these equations have to be solved numerically, but 
simple analytical solutions can be found in two limits. In the long wavelength limit fc <C one obtains 
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2 2 I ^ I 2 



while in the opposite situation, k ^ ujp, one gets 



LOT 



(21) 
(22) 

(23) 



which correspond to the dispersion law of free propagators. Anyway, in the ultraviolet regime the gauge field prop- 
agators receive quantum corrections which are not accounted for by the HTL or kinetic theory approximation. For 
those scales these two approaches arc clearly not valid. 



5 



Now we consider the collective modes which are present in a system described by the tsunami distribution function 
^ alone. These have been considered in Ref. [l^. Upon substituting the tsunami- like distribution function Eq. 
in Eq. ([9]) one deduces that 

<m ^ -4. + ^^^^ - ^f-^^ . m 

^ \ w — k • V (oj — k • v)'' / 

2 _ 

where the plasma frequency of the jet is given by lu?^^^ ~ 

Upon substituting this expression of the polarization tensor in Eq. (|18p . but using the dielectric tensor associated 
with the jet ej^j, one obtains in a straightforward way the dispersion laws of the collective modes of the system. These 

will in general depend on |k|, |v| and cos 6* = k • v. For any value of these parameters we find that the system is 
stable. One can obtain analytical expression of the collective modes for k || v, or k _L v, or w = 1. In the last case we 
find the following non-trivial solutions 

= ojf^^ + fc^ , w = fc cos 6* . (25) 

IV. COLLECTIVE MODES IN THE SYSTEM COMPOSED BY THE QGP AND JET 

We now consider the system composed by an equilibrated QGP traversed by a jet of particles. We are interested 
in very short time scales when the Vlasov approximation can be employed. The effect of the beam of particles is 
to induce a color current, which provides a contribution to the polarization tensor. Therefore apart from the HTL 
polarization tensor, we have to include the polarization tensor due to the jet. In this very short time regime, the 
polarization tensor of the whole system is additive, meaning that 

iir{k) = u>;^k)+u^:,ik) . (26) 

The total dielectric tensor is now obtained considering both the effects of the plasma and the jet 

el^iu;,k)^6^^ + %, (27) 
and the dispersion laws of the collective modes of the whole system can be determined solving the equation 

det \k^S'^ - k'k^ - uj'^el' (fc)! = . (28) 



The solutions of this equation depend on |k|, |v|, cosO = k • v, and on 

,2 



ul^ul+ul, (29) 



and 



(30) 

We will analyze separately the cases where k || v, k _L v and the case of arbitrary values of the angle between k and 

V. 

As we discussed in the previous section, when the plasma and the jet do not interact, they have stable collective 
modes. However, once we consider the composed system of plasma and jet interacting via mean gauge field interactions, 
unstable gauge modes may appear. In the following subsections we will analyze under which circumstances the system 
becomes unstable. 

A. k parallel to v 

Here we consider the case where the momentum of the collective mode is collinear with the velocity of the jet. The 
dispersion law of the unstable mode can be obtained from Eq. (|28p upon setting 9 ~ Q and corresponds to one of the 
solutions of the equation 

iol,{l-v^)-e^{u~kvf ^Q. (31) 
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This equation admits a solution of the form lu = a + iT, corresponding to an unstable propagating mode. We find 
that such an instability is present for any value of the velocity ?; < 1 and for any value of the momentum k. 

In Fig. [T] we present the plot of the imaginary part of the dispersion law of the unstable collective mode as a function 
of the momentum for four different values of the velocity of the jet and for two different values of the parameter b. The 
left panel corresponds to b ~ 0.1 and the right panel corresponds to = 0.02. In both cases the largest value of T is 
obtained for velocities of the jet v ^ 0.9. The corresponding value of the momentum is fc « 2L0t and is approximately 
independent of b. In agreement with the finding of Ref. in the v ^ 1 limit, which corresponds to considering the 
massless limit for the particles of the jet, this mode becomes stable. This can be easily understood from Eq. (|3ip . 
Indeed for v = 1, the dispersion law of the longitudinal modes are not affected by the presence of the jet, because the 
first term on the left hand side of Eq. (|31|) vanishes. Therefore in this limit the jet does not have any effect on these 
modes and the system becomes stable. 

It is interesting to notice that with increasing values of b the value of F at the maximum does in general increase. 
The qualitative reason is that with increasing values of b the plasma frequency of the jet increases and therefore the 
destabilizing effect of the jet becomes larger. Indeed it is clear that setting 6 = is equivalent to neglect the presence 
of the jet at all, as can be easily deduced from Eq. (j3ip . With increasing b (at fixed ujp), the first term on the left 
hand side of Eq. (|3ip becomes larger and the corresponding influence on the solution of the equations grows as well. 
The value of the momentum corresponding to the maximum of F is instead approximately independent of b. 
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FIG. 1: (color online) Imaginary part of the dispersion law of the unstable longitudinal mode for the system composed by a 
plasma and a jet in the case k || v as a function of the momentum of the mode at 6 = 0.1 (left) and at 6 = 0.02 (right) for four 
different values of the velocity of the jet, |v|. 



B. 



k orthogonal to v 



The dispersion law of the unstable mode can be obtained from Eq. (j28p upon setting 9 ^ tt/2 and corresponds to 
one of the solutions of the equation 



-jot 



4 

^jct — T" 



= 0. 



(32) 



We notice that the solutions of Eq. p2p are non propagating, in the sense that they are pure imaginary of the form 
uj = iT. 

We have numerically solved Eq. and reported in Figl^] the plot of the imaginary part of the dispersion law of 
the unstable collective mode as a function of the momentum. 

As for the longitudinal case we find that with increasing values of b, the values of F at the maximum increase and 
the corresponding value of the momentum remains approximately the same. Differently from the longitudinal case, 
there exists a maximum value of the momentum, fcmax, such that for values of the momentum k > /cmax there is no 
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FIG. 2: (color online) Imaginary part of the dispersion law of the unstable mode for the system composed by a plasma and a 
jet in the case k ± v as a function of the momentum of the mode at 6 = 0.1 (left) and at 6 = 0.02 (right) for four different 
values of the velocity of the jet, |v|. 



unstable solution. This quantity can be obtained by solving Eq. ((32|) in the limit F ^ and is given by the expression 



3w2 - 2t>26 - 6 



l-i;2 

Requiring that fcmax is real, we find that there is an unstable mode for 

3t;2 



or equivalently for 



V > Wmin(&) = 



3-26 



(33) 



(34) 



(35) 



Therefore there is a threshold value for the velocity of the jet, that depends on b. If the velocity of the jet is smaller 
than Vmin{b), then no unstable mode with momentum orthogonal to v appears. Notice that the minimum value of 
the velocity vanishes for vanishing values of b. 



C. Arbitrary angles 

We shall now consider various values of the angle between the velocity of the plasma and the momentum of the 
collective mode. The dispersion laws are obtained solving Eq. ([^5]) for arbitrary values of 9. Equation can be 
rewritten as follows: 



(36) 



where 



.4^^ = 1- 



(cj2 - sin^ 

{uj — kv cos 6)'^ ' 
kvsm6 {uj'^ — k'^)v'^ sinO cost 



uj — kvcos.9 (ti> — few cos 0)2 

^2(1-^2 COS2 0) 



(37) 



(w — kv cos I 
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FIG. 3: (color online) Imaginary part of the dispersion law for the unstable mode as a function of k for three different values 
of the velocity of the jet |v| and four different angles between k and v for h — 0.1. The left panel corresponds to v = 0.4. The 
central panel corresponds to « = 0.6 and the right panel to « = 0.9. 



We have numerically solved Eg. pop and we report the results of our analysis in Fig. [31 Here we present the plot of 
the imaginary part of the dispersion law F as a function of the momentum k for three different values of the velocity 
and four different values of the angle between k and v. The plots are obtained for b ~ 0.1. 

For relatively small values of the velocity, the representative plot is shown on the left panel of Fig. [3] for w = 0.4, 
there is no favored unstable direction of the momentum. Indeed the unstable modes corresponding to different angles 
have roughly the same maximum value of F. With increasing values of the velocity of the jet, modes that are not 
collinear with the velocity of the jet are favored. In the central and right panel of Fig. [3] we show the results obtained 
with V = 0.6 and v = 0.9, showing that the largest value of the maximum of F corresponds to the mode with 6 ~ 7r/2. 



V. FLUID VERSUS KINETIC THEORY APPROACHES 



In our previous paper [5| we discussed the appearance of jet-induced plasma instabilities using the chromohydro- 
dynamical approach developed in Q. Here we compare the two formalisms, as well as the results one gets from the 
two approaches for the dispersion laws of the plasma waves in different situations. 

In the Vlasov approximation valid for sufficiently short time scales when the collision terms can be neglected, the 
transport equations can be re- written in a fluid-like form, as first realized for non-relativistic electromagnetic plasmas 
0. In Ref. the same strategy was followed to study the transport equations obeyed by colored particles. Let us 
recall the basic inputs of such a methodology. 

After defining the following two momenta moments of the quark distribution function (antiquarks and gluons are 
treated similarly) 

n^^ix) = I p^Q{p,x) , T'-(a;) = / p'Vg(p,x) , (38) 

one can obtain their dynamical evolution employing the kinetic equations. In particular, after integrating over 
momenta Eq. ((T|) and Eq. ([T]) times a four momentum (with C = 0) one gets 

D^n^^O, i?^T^''-|{F/,n^} = 0, (39) 

respectively. Note that in this short time approximation there are more apparent conservation laws than in a real 
hydrodynamical configuration. For example, Eqs. (j39p suggest that the quark particle density is a conserved quantity, 
while only baryon number (that is, the quark minus antiquark particle density) is conserved at long time scales. 

The fluid approach, originally developed to study the short time evolution of a system, might be used at longer 
time scales if the effect of the collision terms in the kinetic equations could be encoded into new terms in the fluid 
equations proportional to some transport coefficients 0. We will not consider such a possibility here. 

In the fluid approach one considers a specific form of the momenta moments defined above. Assuming that the 
unperturbed quantities take the ideal fiuid form, their small deviations are then expressed as 



= nu'^ + n (5u^ -|- Sn u'^ 



(40) 
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T^"' = (e + p)u'' u^-p g^"" + {5e + 5p)u^' u" + [l + p) (u^ Su'' + Su^" u") ~ Sp g^"" . (41) 

The relevant question now is whether the fluid approach, which essentiafly only studies the dynamical evolution 
of two momenta moments of the distribution function, is describing to an acceptable level of accuracy the physical 
situation one is interested in. The answer to this question depends on the situation one wants to describe, and 
more specifically, on the form of the unperturbed distribution function. In particular we are interested in seeing at 
which level of accuracy the fluid approach reproduces the behavior of the kinetic equations obtained in the HTL 
approximation. We compare in the subsequent subsections the dispersion laws obtained in the two approaches for 
different configurations. 

A. Collective modes for stable configurations 

Let us consider the physics of the beam of particles. In this case we considered that the distribution function is of 
the delta- like form given in Eq. ([8]), that in the non-relativistic plasma literature is known as the "cold beam approx- 
imation" . For relativistic systems this tsunami-like approximation is also correct for ultrarelativistic velocities of the 
jet 0- For the unperturbed distribution function, specifying only the two momenta moments correctly characterizes 
the physics of the system. For their fluctuations it was found in Rcf. that the polarization tensor associated to the 
distribution function ([H]) agrees exactly with the polarization tensor one finds with the chromohydrodynamical equa- 
tions by neglecting the effect of the pressure gradients. One then concludes that in this situation the fluid methods 
describe the color polarization effects in the same way as kinetic theory. 

Let us now consider the equilibrated QGP without the jet. The distribution functions of the particles are those 
of thermal equilibrium, and thus it is obvious that the unperturbed hydrodynamical variables are enough to specify 
the state of the system. Can one describe the color polarization by considering only the first and second momenta 
moments of 5Q1 To answer this question it is instructive to compare the dispersion laws of both the transverse 
and longitudinal collective modes with those obtained in the HTL approximation. Employing the hydrodynamical 
equations of Ref. Q we find that the transverse and longitudinal dispersion laws are given by 

= ujl + k^ , (42) 
u^l = col + ic:fk\ (43) 

where c° is defined by the relation 6p°' = (c°)^(5e° [H. In a fluid approach a relation between the energy and pressure 
variables is needed in order to close the system of equations. In Refs. @ and [H| it was chosen as = 1/ -\/3, assuming 
the conformal limit. However in a fluid description of the plasma one can consider as a free parameter. 

From the behavior of these dispersion laws, one can infer that the hydrodynamical approach is valid for long 
wavelengths, i.e. for k ujp. Indeed in this case the dispersion laws in Eqs. (|42p and (j43p approximately reproduce 
the corresponding HTL dispersion laws. Also the limit fc 3> seems to be approximately well described in the fluid 
approach because the dispersion laws of both the fluid and the kinetic theory approach tend to those corresponding 
to free propagators. However, for values of order A; ~ Wp, we can expect some discrepancies between the fluid and 
kinetic approaches. 

B. Collective modes for the unstable configurations 

Let us now compare the dynamics of the equilibrated plasma traversed by a ultrarelativistic jet of particles in 
the fluid and kinetic approaches. In the two approaches unstable collective gauge modes appear. We will first deal 
with transverse modes that are the dominant ones in the kinetic theory approach. Based on the considerations of 
the previous subsection, we may expect that for very short time phenomena the fluid approach provides the correct 
dispersion laws in both the situations k ^ ujp and k ^ ujp. Wc have made a comparison of the unstable growth rates 
obtained within the two formalisms, and this is indeed what we found. The results of our analysis for the case k _L v 
are reported in Fig. [4] The left panel of Fig. |4] corresponds to = l/\/3, dashed lines (red online) are obtained 
with the fluid approach, full lines (black) are obtained with the kinetic theory approach. We have reported the results 
obtained with two different values of the velocity of the jet; the upper lines correspond to v = 1 and the lower lines 
correspond to v = 0.9. The agreement between the two approaches is quite remarkable. However it is possible to 
improve the agreement employing a different value of c°. In the right panel of Fig. 2] we compare the results of the 
two methods for the same values of the velocity of the jet, but with c" = 1/2 [2l[. The difference between the results 
obtained with the two methods is now of ~ 10% at most. 
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FIG. 4: (color online) Comparison of the imaginary part of the dispersion law of the unstable mode for the system composed 
by a plasma and a jet in the case k _L v as a function of the momentum of the mode at b = 0.1 for v — 0.9 (lower curves) and 
V = 1.0 (upper curves) in the conformal limit (left panel) and for cj = 1/2 (right panel). Dashed (red online) lines correspond 
to the results obtained with the fluid approach; full (blue online) lines correspond to results obtained with kinetic theory. 



Regarding the longitudinal modes, there is agreement between the fluid approach and kinetic theory in the fc ^ 
domain. However, beyond this region we find that the fluid approach largely overestimates the growing rate and 
predicts a threshold value for the momentum. The reason for such a disagreement is due to the fact that longitudinal 
modes are propagating and the solution of Eq. ()31|1 consists of a real and an imaginary component. Therefore one 
cannot reproduce properly such results with one single real parameter c° in the regime A: ~ Wp. However, it might 
be possible to match also this regime if one allows for a non-local relation between the colored pressure and energy 
density. This possibility is discussed in Appendix A. 

There are also two qualitative differences between the results obtained within the kinetic theory method and the 
fluid approach. First of all we notice that in the fluid approach the instabilities develop for velocities w > c° whereas in 
kinetic theory the threshold value of the velocity for the development of the instability is not related to the parameter 
c". This is clearly due to the fact that the equation of state of matter does not enter the kinetic theory picture and 
therefore the colored speed of sound does not play any role. 

Second, in the Vlasov approximation for sufficiently small velocities, v < 0.6, there is no preferred unstable direction, 
whereas for larger velocities the most unstable modes correspond to large angles between k and v. This has to be 
contrasted with the results obtained using fluid equations, where one finds that for velocities v > the most unstable 
mode correspond to momenta k collincar with the velocity of the jet, whereas for ultrarelativistic velocities v < 1 the 
unstable modes corresponding to angles 6* > tt/S are dominant and the most unstable mode corresponds to 6* ~ 7r/4. 



We have studied a system composed by an equilibrated quark-gluon plasma traversed by an energetic jet of particles 
using transport theory. We have assumed that the interaction between the jet and the plasma is only mediated by mean 
gauge field. Moreover the plasma has been considered to be in thermal equilibrium, whereas the initial configuration 
of the jet has been described employing a tsunami-like distribution function. At very short time scales after the 
interaction between the jet and the plasma sets in we find that the effect of the jet is to destabilize the plasma 
producing an exponential growth of collective gauge fields. For relatively small values of the velocity, v < 0.4, we do 
not find that any particular direction in momentum space is favored. With increasing values of v the most unstable 
mode is characterized by having the momentum orthogonal to the direction of the velocity of the jet. 

We have also compared the results obtained in the present analysis with those obtained in Q where fluid equations 
have been employed in the description of the problem. The two formalisms describe qualitatively the same effect: 
the appearance of jet-induced plasma instabilities. At very long wavelength the growth rates obtained from the two 
approaches of the unstable modes match, while they differ for moderate wavelength values. However, the discrepancies 
in this region can be reduced by a proper choice of the equation of state used in the chromohydrodynamical equations. 



VI. DISCUSSION 
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The description of the plasma and of the jet that we have carried out is simplistic in many ways. First of all we have 
assumed that the jet of particles can be described employing a delta like distribution function. It is of course possible 
to employ more realistic distribution function for the particles of the jet, but this would lead to a more complicate 
numerical analysis. We have also neglected possible saturation mechanisms of the gauge instabilities. Unfortunately 
the saturation time scale cannot be evaluated employing a simple linear analysis, but requires numerical studies. 
Moreover at sufficiently long times the effect of collisions should be taken into account, through the addition of the 
corresponding terms in the transport equations. Indeed if the instabilities do not saturate due to non-linear non- 
Abclian effects, the collisions will probably anyway stop their growth. Further, the scatterings of the particles in the 
jet with those of the plasma will provide coUisional and radiative jet energy loss, which however will be relevant at 
time scales longer than the mean free path time. The consideration of all the above mentioned effects is required for 
a complete analysis of the possible contribution of the proposed mechanism to the jet quenching of energetic partons. 

Even with the simplification employed in this paper we believe that the effect we discuss is relevant because 
it provides a new jet energy loss mechanism, based on pure collective effects, which might be important for the 
description of jet quenching in heavy ion collisions. 



APPENDIX A: IMPROVED CHROMOHYDRODYNAMICAL FORMULATION 



Here we discuss how one can improve the agreement between the results obtained with kinetic theory in the 
HTL approximation and the fluid approach of Ref. Suppose that we want to exactly reproduce the behavior 
of the longitudinal modes for the plasma without a jet in the two formalisms. In both kinetic theory and the 
chromohydrodynamical approach the dispersion law of the longitudinal mode is obtained by equating to zero the 
longitudinal component of the polarization tensor. In order to obtain the same dispersion law for the longitudinal 
modes one can define an "effective" speed of sound (we refer to the parameter which relates the colored variables Sp"' 
and ^e") in such a way that the longitudinal component of the dielectric tensor in the HTL approximation and in the 
fluid approach are equal. It follows that the "effective" speed of sound satisfles the equation 

<{y)^J — ^ — + 4, (Ai) 

y 3(1 + ^ log 1-^) y^' 

where y — k/oj. The plot of this function is reported in FigIS] for real values of y, left panel, and pure imaginary 
values of y, right panel. 

li ' ^ ' r , . . . . . . . . . , 
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FIG. 5: Colored Effective speed of sound defined in Eq. (|A1|I . In the left panel real values of the ratio k/u} have been considered. 
The right panel corresponds to pure imaginary values of k/uj. 



For real values of y the effective speed of sound is an increasing function of y and assumes real values that are 
larger or equal -y/S/S (and therefore larger than the conformal value) for 0<?/<l. At uj = k this parameter reaches 
the largest possible value c° = 1. For values of uj larger than k the longitudinal dielectric tensor becomes imaginary 
and therefore it is not possible to reproduce its behavior with a real function. Therefore the disagreement between 
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the results obtained with the two methods in the description of the longitu dina l modes are due to the fact that for 
k <^ uj the value that one should employ is not the conformal value c° = -\/l/3 but c° ^ -^/S/S, whereas for k > uj 
one needs an imaginary speed of sound. 

It is quite interesting that for pure imaginary values of y, corresponding to imaginary values of k or to imaginary 
values of uj, the effective speed of sound is always real. It is a decreasing function of y and assumes values between 
■\/3/5 and -^/l/S, i.e. between the HTL value (in the long wavelength limit) and the conformal limit. Since transverse 
modes are not propagating, they correspond to pure imaginary frequencies and it is then clear that a real speed of 
sound suffices in approximating the behavior of these modes. However it docs not explain why the value that gives a 
good agreement between the two approaches for the transverse mode is c° = 1/2, that has to be determined with a 
different method [2l[. 

It is worth remarking that employing the effective speed of sound defined above means that we are not considering 
the conformal limit and we match the only free parameter of the fluid approach by means of the weak coupling theory. 
This completely sets out the fluid equations. 

We can now compute the dispersion laws of the longitudinal and transverse modes for the system composed by the 
plasma traversed by the jet employing the effective speed of sound deflned in Eq. (jAip . Since the effective speed of 
sound is deflned in such a way that the longitudinal component of the polarization tensor of the plasma in the HTL 
and fluid approximations are the same, it follows from Eq. (|3ip that the longitudinal modes for the system composed 
by the plasma and the jet are exactly the same. Regarding the transverse modes, the results of the fluid approach 
are compared with those of the HTL approximation in Fig [B] 




FIG. 6: (color online) Comparison of the imaginary part of the dispersion law of the unstable mode for the system composed 
by a plasma and a jet in the case k _L v as a function of the momentum of the mode at b = 0.1 for v = 0.9 (lower curves) and 
V = 1.0 (upper curves) employing the effective speed of sound defined in Eq. (|A1|I . Dashed (red online) lines correspond to 
results obtained with the fluid approach; full (black) lines correspond to results obtained with the kinetic theory approach. 



The agreement is quite satisfactory, but is not as good as in Fig. |4l The reason being that the "preferred" value 
of the speed of sound for the transverse mode is c° = 1/2, whereas we arc employing the "effective" speed of sound 
reported in Eq. (jAip . As can be seen from the right hand side of Fig. [5] the values of such a function are between 
y/S/b and \/l/3 and are therefore larger than the value = 1/2. 
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